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I . INTRODUCTION 


The present report describes a new device in which unusual discharge 
characteristics with generation of neutrons from nuclear fusion reactions 
are observed at very modest input power levels of the order of 100 W. 

In addition. X-rays of energy higher than what the 1000 V amplitude of 
the applied high frequency (HF) voltage would normally allow, are generated. 

The two main advantages of a HF discharge in general are: 1) the 

principal possibility of resonance phenomena which may accelerate particles 
or may excite collective modes up to very high energies with modest HF power 
inputs; 2) the absence of electrodes in contact with the plasma. The 
spectacular results obtained with the new device are at least in part due 
to these two major advantages. 

There is, however, another important, fundamentally new advantage 
which is directly related to the name of the new device: ENRAD. This 

acronym means "Electrostatically Confined Nuclear Radiation Device" and 
was coined by the authors on the basis of physical intuition long before 
the novel confining principle of the device, which justifies its name, 
was understood. The new advantage is based on a new confinement principle 
which is stated here for the first time: rather than confine the plasma 
from the outside with a field, it is preferable to confine the field from 
the outside with a plasma, i.e., to have the electric field localized in- 
side a radiation bubble, surrounded by a plasma which in turn is embedded 
in layers of neutral gas. The basic reason is that the field is coherent, 
and, therefore, characterized by fewer degrees of freedom and more stable, 
while the plasma is incoherent and characterized by a larger number of 
degrees of freedom, many of which are very unstable. The plasma does not 
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need to be in thermal equilibrium but must form a resonant state together 
with the field. This combined state of plasma and electrostatic field is 
a new type of solution of a nonlinear equation which describes what is in 
essence a localized standing electrostatic (Langmuir) wave, a particular 
type of soliton, and a special case of what has become known in the last 
decade as plasma caviton. The properly tuned discharge consists of one or 
11101:6 pairs of luminous disks. Each pair is invariably separated by a 
sharply defined dark space of about 0.5 cm width which does not depend on 
pressure and is, therefore, fundamentally different from the familiar 
striations often observed in d.c. or a.c. discharges. The disks occur 
only in pairs, and are the source of a spectrocopically identified axially 
directed beam of deuterons, and also of keV x-rays. 

On the basis of this experimental evidence, we have been able to iden- 
tify the discharge with just one pair of disks as a ground state electro- 
static plasma caviton and the states with n pairs of disks as the corres- 
ponding (n-l)-th excited cavi tonic state. The caviton is resonant with, 
and pumped by, the applied HF field. The sharply defined dark space is 
occupied by the strong localized electric field of the caviton, which 
exceeds the amplitude of the applied HF field by a factor of the order 
of 10 3 and which oscillates at a frequency which differs from the local 
plasma frequency by a small amount determined as the eigenvalue of the 
cavi tonic state. The luminous disks represent the walls of the density 
cavity which houses the localized electric field of the caviton, and are 
active in defining the caviton. Inside the caviton the plasma concentra- 
tion is by about two orders of magnitude lower than the surrounding plasma 
concentration due to the action of the ponderomotive force of the localized 
field which pumps the plasma out and thereby reduces the convective loss 
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of energy of the caviton in form of electrostatic (Langmuir) waves and 
electromagnetic waves. This, in turn, facilitates the further increase 
of the localized resonant electric field and the deeper excavation of the 
density cavity. 

The initial generation of the caviton is explained in detail in the 

present report (Sec. V). The explanation is based on a resonance between 

2 1/2 

the applied HF field of frequency a) and the plasma frequency to = (47rne /m) 

p 

in the center of the device as soon as the HF power is applied to the device 

2 1/2 

and as soon as the electron concentration reaches the level to (m/47re ) in 

o 

the center of the device. The ENRAD cavitons are obtained at pressures of 
the order of 1 torr which is about a thousand times higher than the pressures 
at which cavitons were obtained previously. 

Due to its unique resonance properties, the ENRAD device has the po- 
tential to be developed into a controlled fusion device. As suggested by 
Handel in 1974, -natural plasma cavitons, pumped by huge atomospheric H^O 
masers may already be known as ball lightning. 

The intriguing question which arises is why does a plasma caviton 
arise in this device and not in any other simple d.c. or a.c. gas discharge? 

To understand why, we note that due to the ring-shape of the electrodes 
and their positioning outside the discharge tube, an extended quasi-homo- 
geneous high-electric field volume is present in the center of the device. 

This is an electrostatic analog of the familiar Helmholtz coils. Due to 
the remoteness of all solid surfaces which could produce increased plasma 
recombination, the plasma density exhibits an extended, flat maximum in 
the center of the device, with a quasi-homogeneous concentration of 
electrons and ions which is steadily increasing from zero after the HF 
field is turned on. After 10-20 periods of the HF-fields at the nominal 
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field amplitude, the local electron plasma frequency (4 Tine /m) in the 

center of the device approaches the frequency of the applied field. 

Since the dielectric constant of the plasma (derived in Sec. V) is £ = 1 - 
2 2 

Ujp /w , it approaches zero in the median plane, i.e., at the center of the 
device. On the other hand S = £E = const, in the space between the elec- 
trodes. Therefore, in the resonant volume at the center of the device 
E ->■ oo when n increases to the point where (jo = 03^. Moreover, this strongly 
enhanced electric field amplitude is in resonance with the free electrostatic 
oscillations of the electrons. Indeed, under the influence of the strong 
electric field present only in the median plane of the device electrons 
are depleted on the +E side and accumulated on the -E side. Electro- 
static repulsion, which is the restoring force of the Langmuir waves, causes 

the electrons to rebound after T where T = 2V W • However, since w = 

p/2 P P o 

0) in the center plane of the device, the electric HF field will have changed 
P 

its direction also meanwhile. Therefore, the electrons will be accelerated 
again, and rebound even stronger, like a harmonic oscillator, whose period 
does not depend on the amplitude; this is true for electrostatic waves. 

The electrons are, therefore, continuously accelerated by the electric field 
which is always in phase. The resonance process enhances the average kinetic 
energy and the kinetic pressure of the electrons in the resonance volume 
dramatically, thereby creating a density cavity. The amplitude of the 
electric field and of the electronic oscillation would increase indefinitely 
if it were not bounded by two processes which carry away approximately 
equal amounts of energy: a) the generation of high amplitude Langmuir waves 
which propagate away from the central plane in both directions towards the 
hollow electrode-planes, and b) the acceleration of ions up to energies of 
tens or hundreds of keV, leading to fusion reactions at very modest HF 


powers 
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The creation of the initial density cavity reduces the loss of electro- 
magnetic and Langmuir wave energy due to the formation of a resonator-like 
structure with conducting cavity walls. This reduction in loss causes 
another, additional electric field increase in the caviton and in turn 
a deeper excavation of the density caviton until a dynamic equilibrium 
structure emerges, with very high values of the electrostatic and plasma 
kinetic energy, with vacuum inside, limited by nonlinear and loss phenomena, 
the ENRAD caviton. This is a dark space between two luminous disks. 

The inner side of the disks appears blueish due to the action of primary 
electrons and ions which have higher energies, while the outside regions 
of the disks appear reddish in color due to the lower excitation energies 
of the bombarding particles and of the radiating species in the tails of 
the plasma distribution of the caviton. 

We can also understand the initial large accumulation of electromagnetic 

field energy at the resonant region in the center of the device by follow- 
. “V c ->■ 

ing the field lines of the Poynting vector S = — E x B which represents 
the energy current density. These enter the discharge tube radially, 

~b 

i.e., laterally through the glass walls, perpendicular to E and B. The 
energy flows in with the group velocity v^ = dto^/dk = c/e which is derived 
from the dispersion relation (Eq. 4.6) 


2 2 2 2 
03 = 03 + C k 

o p 


and is therefore v =»£/( 

g 


2 2 

1-0) /0) ) . At the center of the device 03 = 03 

p o p o 


and v^ = 0, i.e., the buck stops there... Energy, again would pile up 
indefinitely because in the electromagnetic energy current the product of 
energy density w and group velocity must be constant in stationary con- 
ditions, as required by the continuity equation 
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HW -*• 

^+V-(wv) =0 
ot g 

i.e. 

-> — >- ~y 

S = w v = const, along the field lines of S. 

g 

Since v •> 0 we have again w + <». This line of thought is equivalent to 

g 

the initial D = £E argument, used above. 

The most amazing property of the ENRAD device is the excellent 
stability of the obtained cavitons against instabilities and collapse. 

This allows the unlimited c.w. operation of the device at the relatively 
high pressures of more than 1 torr which exceeded by at least three orders 
of magnitude any previous observations of laboratory cavitons. Concomi- 
tantly, the electric field and the particle energies are also by about 
3-4 orders of magnitude larger in our device, therefore allowing for 
fusion reactions and x-ray emission. 

We are sure that the pressure and the energy density can be further 
enhanced in the ENRAD device and we believe that this is the way to solve 


the fusion problem. 



II. EXPERIMENTAL SETUP 


1. The ENRAD Principle 

ENRAD stands for electrostatically confined nuclear radiation device. 
Basically, the intention is to contain a field inside a plasma, rather 
than going the usual way where a plasma is confined within some sort of 
field configuration. For the field to exist within the plasma, it has to 
displace the major part of the plasma within a certain region (dig itself 
into the plasma) . The result should be a low density region surrounded 
by a steep density gradient. The resulting plasma cavity (caviton) gives 
charged particles trapped inside an opportunity to obtain energy in the 
applied radio frequency field (RF) until they become too energetic to follow 
the vibrations of the RF field and are ejected from the cavity region. They 
collide with the plasma wall giving rise to excitation, ionization and possi- 
bly some fusion reactions. 

There are probably more than one electrode configuration which will 
shape a field in a way that it can "dig in." The present experiments 
used two cooper rings, placed outside of the discharge tube, so that the 
coupling is capacitive. Electrode distance is about twice the electrode 
diameter. The wavelength of the applied RF is large compared to the elec- 
trode distance or diameter, so that initially a fairly homogeneous field 
can be expected between the electrodes. 

In the following two different discharges which are based on the ENRAD 
principle are discussed. Both use the same electrode configuration. One 
discharge is a pulsed DC- high voltage discharge, while the other one is 
a low voltage RF discharge. 

The work concerning the high voltage discharge was done more than 
10 years ago, while the work on the RF discharge was done under the present 
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effort. However, only now, after having the results of the RF discharge, 
can we understand the results of the older work. 

2 . High Voltage Discharge 

The high voltage discharge is excited by a voltage pulse, 5 ysec wide 
and 60 kV high. Its repetition rate is 50 hz. The voltage pulse is 
capacitively coupled to the plasma, using the ENRAD configuration. The 
relatively high repetition rate has several advantages. One of them is 
the outgassing of the structure which results from repeated electrical 
pulses, driving out contaminating wall gases and permitting the operation 
of a stable relatively reproducible plasma. Another advantage is that 
the residual gas pressure can be stabilized after a short time of operation. 

The pulse generator used can produce voltage pulses as high as 100 
KV and currents as large as 100 A. The pulse duration is approximately 
five microseconds and the repetition rate may be either 60 Hz or 400 Hz. 

The shape of the pulse is, to a large extent, determined by the load char- 
acteristics which in the case of the gaseous discharge vary considerably 
with the operating condition. Since the high field gradients produce a 
large amount of ozone in the surrounding atmosphere, the tube and elec- 
trodes have been enclosed by a large steel tank equipped with suitable 
ports for observation. This tank serves to confine the ozone and prevent 
its entrance into the experimental environment. It also acts as a good 
X-ray shield. Still another advantage of this enclosure is the electrical 
shielding which results, reducing interference to measuring equipment 
in the immediate vicinity of the discharge. 

Figure II-l shows the schematic diagram of the power supply. The 
trigger circuits consist of an adjustable phase, synchronous pulse 
generator which serves to trigger an intermediate level trigger, a thyratron 



Figure II-l Schematic Diagram of Power Supply and Discharge Tube* 
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circuit which, in turn, triggers a high-power ignitron circuit supply- 
ing energy to the high-level pulse transformer which supplies the 100 KV 
pulse. 

Energy for the high-level pulse is obtained by discharging a trans- 
mission line which can be adjusted to produce various wave forms of the 
resultant pulse. This energy storing line is charged by a diode and 
fired during the nonconducting portion of the cycle thus establishing 
synchronous operation. The pulse repetition rate is thus established by 
the source frequency. Such modulators are employed in most high-powered 
pulse work. 

A switch is provided to change the pulse polarity. In the following 
discussions, the term "positive" or "negative polarity" is defined as 
the potential of the ungrounded electrode. 

The discharge tube consists of a cylindrical Pyrex glass tube 200 
cm long and 8 cm in diameter. Three cylindrical Al-electrodes are placed 
on the outside of the tube as can be seen in Figure II-2. This figure 
also shows the above mentioned enclosure. As pointed out, the power input 
is transferred to the plasma by capacitive coupling. The central elec- 
trode has guard rings in order to reduce the corona discharge on the out- 
side of the tube. The tube was filled with either hydrogen or deuterium 
at pressure levels between 10 and 200 microns. 

A closer look at the relationship between current and voltage reveals 
an important feature of this discharge, namely that the current rings in 
the 1 Mhz range, although the applied voltage, is a dc pulse. This can be 
seen in Figure II- 3 where the applied voltage pulse is +50 KV high and 
4 ysec wide, while the current rings (positive and negative values of 
current) at high frequencies. 
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Figure II-3 Relation between Current (Top Trace) and Voltage (Botton Trace) 
at 50 kV Pos. Pulse Polarity (Time Base 2 psec/cm) for various 
pressures. 
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3. Low Voltage Discharge 
a) Description 

The low voltage system consists of a test section, a matching circuit, 
a driver unit, a gas handling system-pump station and a diagnostic unit as 
shown in Figure II-4. The driver unit is a 500 watt RF power supply (0.5 
to 1.7 Mhz) feeding a matching circuit instead of a conventional antenna 
system. The matching circuit is a variable L-C network which allows match- 
ing of the output impedance of the RF power supply to the imput impedance 
of the test section. 

The test section consists of a 80 cm long glass tube (7740 borosil- 
licate glass) having an outer diameter of 10 cm. According to the ENRAD 
Principle, the HF is capacitively coupled into the fill gas of the test 
section by two ring electrodes which are again attached to the outside wall 
of the glass tube. Intra-electrode distance is typically 5-50 cm. 

The gas fill for the test section is introduced through the gas 
handling system while the vacuum is monitored and maintained through a 
separate port. This configuration allows studies with either static or 
dynamic gas fills. 

b) Operational Characteristics 

Like all gas discharges, this device has a breakdown voltage which 
is a function of pressure and electrode separation. It is usual to pre- 
sent such relationships in the form of E/p vs pd. (E: electric field, 
p: pressure, d: electrode distance) . The breakdown E/p is shown in Figure 
II-5 for a frequency of 1.34 Mhz, which is the primary frequency used for 
the data in this report. The dependence on frequency is weak. As a 
comparison, an argon HF-discharge at 3000 Mhz is shown. [II-l] The dis- 
charge characteristics are quite similar despite the large difference in 
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Fig. II. 4. The ENRAD System 
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Fiq. II-5 . 


E/p vs pd for ENRAD at 1.34 Mhz 




Fig. II-6 . Typical pair of plasmoid disks in Hydrogen. 
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Fig. II-7. 


Two Pairs of Disks in Hydrogen. 



Fig. II-8. Three Pair of Disks in Hydrogen. 





Fig. II-9. One Pair of Disks in Helium. 
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Fig. 11-10 


Number of disk pairs vs pressure for ENRAD device 
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frequency, especially at larger pd values. It seems that after the cycle 
time is shorter than the effective recombination time or diffusion time 
(whichever is shorter), the plasma is stationary and an increase in the 
frequency has very little effect as far as maintenance of the discharge 
is concerned. (This is not to say that other properties such as neutron 
emission are independent of frequency. ) 

Once the discharge is ignited, particular E/p ratios will establish 
certain well defined visible plasma structures. These structures are 
pairs of disks in the following referrred to as "cavitons" which are 
separated by a very well defined dark space. (The separation appears 
black to the eye. ) A pair of disks (a caviton) will exist over a certain 
pressure range (the applied voltage is kept constant) and if a certain 
critical pressure is exceeded, the plasma transforms (abruptly) into 
another mode consisting of 2n disks, or n cavitons where n is an integer. 
Again, each pair is separated by a very well defined dark space. 

Figure I 1-6 shows such a caviton (n=l) while Figure II-7 and II-8 
show two and three cavitons in deuterium. If helium is used, cavitons 
still can be obtained but are not as well defined as in deuterium or 
hydrogen. Figure II-9 shows such a caviton in helium. In any case the 
location of the cavitons is always symmetrical to the electrodes. 

The cavitons should not be mistaken with the familiar striations 
(moving or stationary) of d-c or a-c discharges. The difference is that 
in these discharges the distance between striations changes when the 
pressure is changed; while for the cavitons described here the distances 
remain constant even if the pressure is changed until the critical pressure 
is reached where abruptly the next higher number of disk pairs is formed, 
having then different distances from each other. This relationship be- 
tween pressure and number of disks can be seen in Figure 11-10. 
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If the pressure is too low for formation of cavitons another phenome- 
non can be observed, which is referrred to as degenerated cavitons. 

The cavitons are no longer suspended between the electrodes but are 
now attached to them. It appears that the one wall of the caviton is 
now a crescent, while the other wall is formed by the sharp edge of the 
electrode. Figure 11-11 shows a photograph of such a degenerated caviton. 

Since there is a minimum pressure for formation of regular cavitons, 
an obvious questions is — if there is also a minimum frequency? There- 
fore, the minimum frequency required to form a disk for a given pd- value 
was determined. In Figure 11-12 this relationship is shown. The elec- 
trode distance (d) was varied while the pressure was kept constant (p = 

0.6 torr) . At any frequency below the straight line in Figure 11-12, 
no cavitons will be formed, but a regular diffused discharge will appear. 
However, operation above the straight line will have little effect on 
the then established cavitons. This is confirmed up to 1.7 Mhz. 

Another important piece of information relates to the question if 
the plasma remains in existence during the time when either voltage or 
current is zero. If the plasma should cease to exist at that time, re- 
start of the discharge would be required each time. In this case, con- 
finement of the plasma would not have been achieved. 

Demonstration of a continuous plasma can be done by determining the 
relationship between light output, and electron density. The light out- 
put indicates the presence of excitation rather than ionization. However, 

—8 

lifetimes of excited states are typically 10 seconds, two orders of 
magnitude shorter than the cycle time. Therefore, if the electron density 
were to be lowered significantly or should totally disappear, light out- 
put would cease between cycles since electrons are undoubtedly responsible 
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Fig. 11-11. 


Degenerated Caviton. 
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Fig. 11-12. Electrode separation vs Plasma structure cut-off frequency 
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for excitation of the atoms* Therefore, we believe that the intensity of 
the observed Balmer lines is a good indicator for the presence of electrons* 
Figure 11-13 shows a series of oscilloscope traces depicting total 
luminosity versus time for different operation conditions. Figure II-13a 
is the degenerated caviton mode. Here the lowest light output is observed. 
However, it is fairly constant with time, seemingly very well confined. 

(The signal is negative in all oscilloscope traces in Figure 11-13.) When 
more cavitons are added, several types of oscillations can be observed. 
However, none of these oscillations have a frequency which is equal to or 
harmonic of the applied frequency. The time deflection on the oscilloscope 
traces is: 5 msec/cm. The conclusion which can be drawn from this figure 

is that the light intensity never equals zero, confirming that the plasma 
is in existence at all times. 


increasing increasing 





Fiq. II 13. Total Luminosity Output of ENRAD Device. 
Time Deflection 5 msec/cm. 





III. EXPERIMENTAL RESULTS 

1. X-Ray Emission 

In both discharges emission of x-rays took place. While this can be 
ejected in the case of the high voltage discharge, it is not immediately 
evident that this would also happen in the case of the low voltage dis- 
charge. 

The x-ray emission from the latter was, indeed, very weak and since, 
at the same time, neutrons were observed in the discharge, we cannot 
state beyond a reasonable doubt that the observed TLD readings and x-ray 
film exposures were actually caused by x-rays and not by a y-background. 

Since information on x-ray emission from the low-voltage discharge 
is so marginal, we will discuss first the x-ray emission from the high 
voltage discharge. Since the geometrical set up is very similar, it 
should give a good indication what kind of x-ray emission can be expected 
from a caviton forming electrodeless RF discharge. 

a) The high voltage discharge is a strong source for x-ray radiation. 
This radiation seems to appear at or after the onset of the high frequency 
oscillations of the current and is, therefore, somewhat delayed in re- 
spect to the voltage pulse. This can be seen in Figure III-l where the 
upper trace is the x-ray radiation and the lower trace is the applied vol- 
tage. This figure can be compared to Figure II- 3 where the current trace 
can be found. Careful comparison shows that the x-ray radiation starts 
after the first full swing of the current oscillations. 

The location of the x-ray emitting areas within the discharge tube 
was determined by wrapping a Kodak- Medical X-Ray Film enclosed in a light 
tight envelope around the discharge tube in the region between the elec- 
trodes. Figure I I I- 2 shows the location where the film was placed, while 
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Figure III-l Radiation and Voltage (common time base 2 ysec/cm) 

Upper trace: Radiation 

Lower trace: Voltage 
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FILM ROLLED AROUND TUBE 
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Figure III 2 Location of X-ray Film on Discharge Tube. 
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a more detailed sketch of the whole device, including its enclosure, can 
be found in an earlier figure (Fig- II-2). The results (positive) of the 
wrap around is shown in Figure III-3. The most interesting features are 
the crescent shaped exposure patterns and the sharp transition from exposed 
to unexposed areas, although the long side bars are also a phenomenon 
which needs to be explained. In Figure III- 3 some of the crescents over- 
lap, since the x-ray film overlapped when it was placed. Therefore, a 
positive of the x-ray negative, namely the one shown in Figure II I- 3, was 
rolled up and glued together. The unexposed part was cut away. Figure 
III-4 then shows a photograph of this composite. As one can see there 
are 3 symmetrically aligned crescents, quite similar in appearance as the 
one shown in Figure 11-11, although the latter figure is a photograph of 
the visible glow of the low voltage discharge. 

The pulsed high voltage discharge has a feature which the low voltage 
RF discharge does not have. Namely, since a dc voltage pulse, which is 
applied, can run the ungrounded electrode either positive or negative. 

As it turns out, this makes a difference as far as x-ray production is 
concerned. To show this the x-ray radiation was observed side on using 
a collimator and scintillation detector. Figure III- 5 shows the general 
outlay of the discharge tube, the shielding enclosure and the location 
of the detector. Scans parallel to the axis of the discharge tube were 
made and it was found that the maximum radiation emissions occur close 
to the negative electrode, which is also the electrode where the crescents 
were observed. Figure III-6 shows the scan with the center electrode 
negative, while Figure III- 7 shows the scan with the center electrode 
positive, which means the grounded electrode is negative. As can be seen, 
the radiation maximum is always toward the negative electrode. Its loca- 
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Fig. I I 1—3. Wrap around x-ray from older ENRAD like device. 






Fig. III-4. Three Dimensional Sketch of X-Ray pattern Around Tube. 
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tion is at a certain distance (d) away from the negative electrode. This 
distance is a function of pressure, as shown in Figure III-8. The higher 
the pressure, the closer it is to the negative electrode. Observations 
of the x-ray emission were also made end-on and typical x-ray positives 
obtained are shown in Figure III-9. The x-ray films were placed on the 
tube x-axis 1, 3 and 5 feet away from the window. The exposure which 
shows an intense inner core was 1 foot away while the ones which show an 
unexposed inner core were further away. If one would assume that the 
x-rays come from a layer on the inner surface of the tube (or from the 
crescents close to the inner surface of the tube) , one can explain Figure 
III-9. This is done in Figure III-10 and Figure III-ll. The distance of 
the x-ray film is plotted here in feet, in axial direction while the size 
of the pattern is plotted in mm in radial direction. In this way one 
can triangulate back to the x-ray emitting layer. The location of this 
layer agrees with the side on scans of Figure III-6 and III-7. Its struc- 
ture, of course, can be seen from the wrap-around in Figure III-3. 

Of course, if an agreement with the side on scans exists then also 
a switch in polarity must show up in the end-on x-rays. This is, indeed, 
the case and the difference can be seen by comparing Figure III-6 to 
I II- 7. Indeed, the x-ray emission is moved closer to the negative elec- 
trode. 

B. X-ray Emission from the Low Voltage H-F Discharge 

Since the x-ray emission from this discharge is considerably reduced 
in intensity, no clear determination of its structure could be made. How- 
ever, from the above described results, one can roughly estimate the shape 
and structure of the expected x-ray exposures. 





Fig. II I- 9. End View X-Rays from Older ENRAD-like Device. 
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Fig. III-9. Cont. Shows a Shadowgraph of a Film Holder. 




Figure III-10 X-ray Geometrie (50 kV pos. pulse polarity at 7 y pressure) 
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Since x-ray films are a good dielectric, it was not possible to wrap 
the discharge tube. The RF would heat the x-ray film and melt it. 
Therefore, a film was placed next to the tube (parallel to the tube axis) 
without bending it about 1 cm away from the surface of the tube. The re- 
sultant esqposure (negative) is shown in Figure III-12 which has to be 
compared to the wrap-around (Figure III- 3) . Since the exposure is very 
faint, only superficial conclusions can be drawn. Figure III-12 was taken 
with a caviton in the center of the discharge tube. There is, indeed, a 
region of increased x-ray emission in the center of Figure III-12. The 
other two maxima at the electrode must correspond to the (smaller) maxima 
on the positive electrode, which appear in Figures III-6 and III-7. 

For the end-on exposures, one would not expect a pattern as in Figure 
I I I- 9, since the discharge tube was all glass and did not have the stain- 
less steel neck which was responsible for the inner core and outer zone 
definition, as seen in Figure III-9. Therefore, one really could only 
expect a diffuse disk at best or just a general fog if the radiation is 
isotrop. In contrast to this expectation, the exposures shown in Figures 
I 11-13 through III-17 show a distinct structure. There is, indeed, a 
central blob which could be caused by a collimated beam in Figure III-13. 
The exposures are — as can be seen — reproducible only to some degree, 
they vary in intensity and location. Also the streaks intersecting the 
crescent in Figure III-15 vary from exposure to exposure. The small 
crescent of Figure III-14 may correspond to the most heavily exposed por- 
tion of the blob of Figure III-13. The fact that there are two crescents 
on Figure III-16 is understandable, since these are long range exposures 
(usually 8 hrs.) and the discharge may have been turned off once and re- 
started or the pressure may have changed. This may cause the beam to 
deviate somewhat and so cause the double exposure. 



Electrode 



Figure III* -12 X-ray from electrode region of E1JRAD device 
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Figure III-13 


End-on X-ray Output from RF Discharge 
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Fig. III-14. End-on X-ray Output from RF Discharge. 
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Fig. 


III-16 


End-on X-ray Output of RF Discharge 
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The fact that the exposures are not in the center of the figure indicate 
that the beam does not come out straight but is (unpredict ably) deviated 
to begin with. 

Whatever the case may be/ if the exposures shown in Figures III-13 
through 111-17 are caused by x-rays, these x-rays must be emitted in a 
preferential direction. 

An attempt to measure the intensity of suspected x-ray emission was 
made. A Victoreen model 440 RF/A survey meter was placed next to one 
of the electrodes. The 440 RF/A is a highly sensitive ionization chamber 
with a broad energy response from 12 keV to 1.2 MeV. It is designed to 
perform x-ray dose rate measurements in high RF-fields. It is suited for 
soft x-rays, however, its response time is under these conditions as 
long as 12 seconds. Reproducible readings in the range of 1 mR/h to 3 
mR/h were obtained. However, the long response time does not allow any 
cross correlation with actual operation parameters. 

Also TLD chips were placed next to the tube. These chips were ex- 
posed for three hours and then counted and compared with control chips 
that had not been exposed. The conclusion was that the chips must have 
been exposed to a radiation field of 3.3 mR/h. 

Also, an attempt was made to observe the x-ray pulses individually. 

For this reason a Nal-scintillation detector was placed close to one 
electrode. The output was displayed on an oscilloscope, and is shown 
in Figure I I 1-18. Since the x-ray pulses come randomly, the camera 
shutter had to remain open for several seconds. During this time the dis- 
charge went through several million cycles. Power lines inside (the RF- 
shielded building) picked up the radiowaves and leaked it into the oscillo- 
scope. However, the discharge is very repeatable. E.g., the traces in 
Figure 11-13 are also second long exposures, which means also several 
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Fig. HI-18. 


Individual X-ray Pulses. 
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million repetitions. Therefore, the traces coming from electrical inter- 
ference can be seen as a band at the lower part of Figure III-18. The 
individual x-ray pulses are seen in the upper part of the picture. 

2. Spectroscopic Studies 

a. Experimental Arrangement 

The system used for study of optical emission coming from the plasma 
is shown in Figure III-19. It consists of a 1 m Czerny-Turner grating 
spectrograph/ spectrometer and an optical imaging system. For spectro- 
meter operation, the scanning of the grating is accomplished by a stepping 
motor. Photoelectric registration is accomplished by an EMI 8575 photo- 
multiplier tube and a picoammeter . The particular Czerny-Turner mount 
used is convertible to photographic registration by rotation of an auxili- 
ary mirror into the spectrograph position. 

The Czerny-Turner mount is equipped with a 1200 line/mm grating 
blazed for 4000& (1st order) . Demonstrated resolving power for this in- 
strument is better than O.lR at 400oR using 1st order. 

The slit width used to achieve this resolution is 16 jim. For reasons 
of speed of observation in the studies described in the following, larger 
slit widths were used for photoelectric registration (up to 50 jim) . 

Slit illumination was accomplished with a 1:1 imaging system, placing 
the image of the plasma into the entrance slit plane allowing spatial 
resolution of the observed plasma structures - 

b. Survey Spectra 

Figure II I- 20a shows a photographic spectrum of the emission of the 
ENRAD device operating in hydrogen. The drive frequency was 1.35 MHz and 
the pressure was 0.6 torr. Since a slit width of 25 ym was used a long 
exposure time was required (6 hrs.). 
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Figure III-19 


System for the study of ENRAD optical emission. 
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Figure II I- 20b shows the spectrum of the emission of the ENRAD device 
operating in deuterium under the same condition. Both spectra were taken 
using side-on observation. They show expected features like the , Hg and 
Hy lines in the hydrogen discharge. D a , Dg and Dy can be seen in the 
deuterium case. In addition, features resulting from the emission of the 
hydrogen molecule are very prominent. 

The Stark broadening of the hydrogen lines is very moderate. Based on 

12 

the observed half widths, the electron density has to be less than 10 

3 

electrons/cm . From this, it has to be concluded that even at the relatively 
low gas pressure under which the discharge is operated, the plasma is only 
partially ionized. This conclusion is corroborated by the presence of sub- 
stantial H^-emission. Usually in fully ionized hydrogen plasmas, molecular 
bands are weak or non-existent. For comparison, a spectrum of a hydrogen 
d-c arc is shown. (Figure III-21) This arc was operating at about 15 000 °K 
which means it was almost fully ionized. There are no molecular bands 
visible and the Stark broadening of hydrogen lines is very prominent. 

As pointed out before in the ENRAD device, electron formation takes 
place in the volume rather than in the vicinity of an electrode as is the 
case at discharges having internal electrodes. These latter discharges 
develop cathode and anode fall regions, which are vital for the maintenance 
of the discharge. In the cathode fall region acceleration of the electrons 
leaving the cathode takes place, requiring large electric fields (the major 
portion of the applied voltage is dropped across the cathode fall region) . 
These fields are partially balanced by an electron cloud adjacent to the 
cathode fall region - the space charge. Here a high electron density and 
consequently considerable Stark broadening can be expected, while in the 
positive column, the electron density is moderate. 


m 
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Figure I.II-20a Hydrogen Spectra from ENRAD device. 



Figure III-20b Deuterium spectra from ENRAD device 


Figure III-21 Hydrogen D.C. 
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Figure III-22 is a spectrum of such a discharge. In the proximity 
of the electrodes, some of the spectrum lines (those which are subject to 
Stark broadening) are severely broadened, while only moderate broadening 
is observed in the plasma columns (the center portion of the lines) . The 
two sample spectra of internal electrode discharges have to be compared 
to spectra obtained from the ENRAD device. The comparison shows that over 
the spatial region covered, there are no electron density gradients, 
especially no fall regions, a fact which was already mentioned as one of the 
advantages of a RF discharge with external electrodes. 

c. Intensity Evaluation 

The information contained in spectral emission rests in intensity, 
wavelength and line width. The latter two were discussed above to some 
extent. 

The intensity of a spectrum line is given by 


I = An* hv 


(III.D 


where A is the transition probability, hv is the photon energy and n* is 
the number of atoms which are excited into the upper level of the line of 
interest. 

If equilibrium exists or is at least locally approximated (local 
thermodynamic equilibrium, LTE) the quantity n* can be computed using the 
Boltzmann factor: 
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X 

kT 


P 

kT 


(III. 2) 


where n Q is the total particle density (p: pressure) g is the statistical 
weight (degeneracy) ; U is the partition function, E ^ is the excitation 
energy, and kT is the thermal energy. 
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Figure HI-22 Sectrum showing Stark broadening in electrode region. 
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If the relationship for the intensity of a spectrum line 


I = A h rn 

A o U 


g - kT 
e 


(III. 3) 


is rewritten in logarithmic form 
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(III. 4) 


Using this relationship, one can relate the observable intensity to the 
excitation energy in a straight-forward way. If the factor where 

the intensity is the only quantity which needs to be measured, is plotted 
on a logarithmic scale versus the excitation energy, a straight line with 
the slope - should be obtained. Such a plot is usually referred to as 
a "Boltzmann Plot." 

Any deviation from the straight line indicates deviation from a Maxwell 
Boltzmann distribution of states (which is the prerequisite for the valid- 
ity of the Boltzmann factor) . A point laying above the line suggests an 
overpopulation of the particular state, while a point laying below the line 
suggests underpopulation of this state. 

A resonance phenomenon is a non- random event. Therefore, one would 
expect that for a resonating laser cavity, at least the point belonging 
to the upper laser level would not fall on the straight line. The same 
has to be true for the ENRAD device if it is truly a resonance machine. 

For this reason, two Boltzmann plots have been taken. One was taken under 
conditions where the ENRAD device was operated like an ordinary RF dis- 
charge, showing no particular plasma structure. The resulting Boltzmann 
plot is shown in Figure III-23. The deviation of the points from the 
3400°K line are probably within the error of measurement. The conclusion 
is that the plasma is more or less in LTE at a moderate temperature of 



Log(lx/gA) 
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Figure III-23 Boltzmann Plot of the Hydrogen H^, Hg, H^, and Lines 
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340Q°K. This temperature will produce only a small thermal ionization 
of the plasma, which is consistent with previous conclusions. 

However, if the ENRAD device is operated in resonance conditions, 
where x-ray and neutron emissions are observed, the Boltzmann plot is 
changed drastically. Figure III-24 shows such a Boltzmann plot for the 
deuterium case. A substantial deviation from LTE can be seen. In order 
to appreciate the magnitude of the deviation, it is convenient to super- 
impose Figure III-23 and Figure III-24. To be able to do this properly, 
one should appreciate that the intensities usually used for Boltzmann 
plots are relative intensities. A closer examination of equation III-4 
shows that only relative intensities are required. Therefore, the inter- 
cept of the straight line with the ordinate (zero excitation energy or 
ground state population) is only correct within a constant factor which 
is different for each gas. Therefore, in order to superimpose both figures, 
one of them has to be moved parallel to the ordinate until one point of 
one figure matches the corresponding point on the other. In case of 
Figure III-25, the point for Dy (Hy) were matched. As can be seen from 
this figure, the deviations are large, covering several orders of magni- 
tude. Although the decision to match the point for Dy(Hy) is arbitrary; 
there is a good argument for why they should be matched in contrast to 
the points H a (D a ) . The slope of a line drawn from Dy to Dg is roughly 
the same as the slope of the Boltzmann Plot for under non-resonant 
conditions. (Figure I I 1-2 3) . It is reasonable to assume that both con- 
ditions produce about the same electron temperatures (for such electrons 
which are thermalized) . Therefore, the conclusion drawn from Figure III-25 
is that D a is heavily underpopulated when the device is operated with 
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Figure III-24 Boltzmann Plot of Deuterium D a> Dg, and D^ Lines 
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cavitons established. If one assumes that the walls of the cavitons 
are formed by ions interacting with the RF field, the observed radiation 
has to come from excited atoms, which are freshly recombined ions. Since 
the cross sections for recombination into higher laying states are larger 
than into lower laying states, it is understandable that the higher laying 



which explains why this level is underpopulated. 


d. Doppler Shifts 

If the ENRAD device indeed produces directed ion velocities — namely 

a beam impinging on the caviton walls — this fact should be detectable by 

measurements using the Doppler effect. However, if the high ion velocity 

is random (temperature), a Doppler broadening of the Balmer lines should 

be detectable. If the beam should exhibit some spread in energy without 

being randomized in direction (or only randomized to a small degree) , a 

broadened Doppler shifted line is expected. 

To demonstrate an existing Doppler shift side-on and end-on light was 

superimposed on the spectroscopic plate, without moving the peak between 

exposures. The optical arrangement for doing this can be seen in Figure 

III-19. A blue shift for the end-on radiation is clearly visible in the 

photographic spectrum depicted in Figure III-26. The size of the shift 

is about 1?, which translates into 40 eV directed motion. 

The end-on observation spectrum shows also a feature at a part of 

the D 0 line which corresponds to the part of the plasma which is in the 
P 

vicinity of the wall. This feature is an additional broadening, which 
is shaded towards the blue. However, the information available does not 
allow a distinction whether this feature is an additional Doppler shift or 


a Stark broadening 
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Combined Boltzmann Plots of Hydrogen and 
Deuterium Showing Magnitude of Inversion. 



Fig„ III— 26. Doppler Shift for and Side On vx End On 
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3. Neutron Emission 


In the following we will show that the ENRAD device emits neutrons 
when operated in deuterium and when cavitons are present. 

The results presented here were obtained under this effort only. 

Since the neutron emission is quite unexpected for such a low voltage dis- 
charge, we felt it was necessary to have the neutron emission verified 
by a researcher not involved in this project. This independent measurement 
was done by Dr. William H. Ellis, University of Florida, and is reported 
in Appendix A. 

Figure I II- 27 shows the system used to study the emission of neutrons 
from the deuterium plasma structures. The detector is a BF^ tube operating 
in the proportional counting mode. The tube and associated electronics 
are RF shielded and extensive tests were performed to guarantee that the 
RF interference did not influence the counter system. The system was then 
calibrated so that an exposure of 3.6 x 10 ± 20% neutrons would result 

in 2000 counts or 1800 ± 350 neutrons per count. The same moderator was 
used for calibration and all tests. 

The system was filled to 300 torr of deuterium and pumped down to 
50 microns within 10 minutes. During the 10 minute cycle, the neutron 
detector was turned on and counts were accumulated. 

Table III-l gives the results of the various runs with the neutron 
counting system. With D^ gas, statistically valid results show that 
neutrons are being produced by the plasma. These neutrons are expected 
to result from the 
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Figure IXI-27 System for the study of neutrons from ENRAD device. 
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TABLE III-l 


Gas 

Counts 

Net 

Counts 

with moderator 

115.20 

17.25 

78.45 

20.03 

w/o moderator 

37.67 

7.23 


- 

with moderator 

42.50 

4.95 

5.75 

15.73** 

w/o moderator 

30.00 

3.87 


- 

*Air with & w/o moderator 

40.5 


- 

*Evacuated Tube with & w/o 
moderator 

33.0 


- 

Background 

36.75 

10.78 




*Average of * items) 

**Not statistically valid 
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Neutron emission from the ENRAD is low. Therefore, to verify that 
the detection system was not influenced by RF interference, an alterna- 
tive technique of neutron detection, namely neutron activation, was per- 
formed. 

A gold coin enclosed in a paraffin moderator was irradiated for 24 
hours. The pd (pressure x interelectrode distance) was such that 6 pair 
of disks were formed. The gold coin was counted immediately after irradi- 
ation. The result was a corrected count rate of 80 CPM. Figure III-28 
shows the decay curve for the foil. The slope of the decay curve computed 
to be a half-life of 2.7 days which is the half-life of gold. This in- 
dicated the BF^ neutron counter was functioning correctly. 

The emission of neutrons by the ENRAD device has several parameters 
which were measured. The first parameter measured was the neutron count 
vs frequency. Figure I 11-29 shows a plot of the data for frequencies be- 
tween 1.20 and 1.45 megahertz. The pressure (with a constant d) varied 
during the taking of the data. It varied from a high value sufficient 
to produce a diffused discharge to a low value capable of producing the 
degenerated disk discharge within 15 minutes. The data does not indicate 
that the shape of the fitted curve needs to be a straight line, but may 
be a resonance shape with a peak about 1.34 megahertz. 

Another parameter measured was the location of the maximum neutron 
emission. The results indicate that the location for maximum neutron 
emission spreads from around the mid-point in the discharge to the grounded 
electrode. The rotational position about the ENRAD axis had a homogeneous 
neutron emission. 

The last parameter measured was the amount of neutron emission vs 
percent nitrogen admixture. It was found that .5% nitrogen admixture to 
the deuterium increased the neutron output without destroying the plasmoid 
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disks. Figure I I I- 30 shows a plot of this data. As the admixture went to 
zero, the neutron count rate stabilized at approximately 300 CPM. 

The admixture altered the stability of the plasmoid disks. An in- 
stability manifested by flickering between n and n + 1 disk pairs was ob- 
served. With an admixture of nitrogen, the flickering between n and n + 1 
disks could be maintained for several seconds. This is not to say that 
during the instabilities the neutron output increased, since the neutron 
counting system does not resolve short times (the order of seconds) . How- 
ever, these results could be an indication that the neutrons are genera- 
ted by instabilities. 
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Fig. III-30. Neutron Count Rate vs Percent Nitrogen Admixture 




XV. EXPERIMENTAL OBSERVATIONS OF CAVITONS 
1. Probe Measurements 

Cavitons have been observed before by others at lower pressures. 

In reference IV-1 the time averaged charge densities around such cavitons 
are reported. These measurements are made with an ion beam probe. Al- 
though the RF discharge used in reference IV-1 had electrodes, while ours 
is electrodeless, certain similarities exist, e.g., they also report 
observation of a dark sheath. Figure IV-1 shows some of their results. 

The striking feature is the charge walls surrounding the cavitons. 

We like to compare these (very sophisticated) measurements to some 
old measurements of our own made with a Langmuir probe (Figure IV- 2 and 
IV-3). The Langmuir probe disturbs the plasma to some extent. Also, the 
interpretation of the probe data assumes a Maxwellian electron distri- 
bution. For these reasons our probe measurements have to be regarded 
as somewhat crude and they are only reported here because of their simi- 
larities with the measurements in reference IV-1. 

-2 

Figure IV-2 shows the results for 1.5 x 10 torr. At this pressure 
only degenerated cavitons (see also Figure 11-11) exist. The electrode 
to which those cavitons attach is located at 12 cm. An indication of a 
surrounding charge wall is visible. The electrons in this wall region 
cannot be expected to have a Maxwellian velocity distribution, since they 
are accelerated towards the wall. The Langmuir probe evaluation process, 
therefore, overestimates their temperature. 

Figure IV-3 shows the results for 6 x 10 ^ torr. The electron density 

shows now a minimum in the center and moderate walls towards the electrode. 

-2 

Results in Figure IV-2 (1.5 x 10 torr) have to be compared with the 

-4 -1 

1.0 x 1. curve in Figure IV-1 and the result in Figure IV-3 (6 x 10 torr) 
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Figure IV-1 Charge Density Walls Abound Cavitons. After Ref. (IV-1). 
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-4 

have to be compared with the 3 x 10 torr curve in Figure IV-1. Ob- 
viously, the Langmuir probe measurements cannot bring out the steep walls 
of the caviton as good as the ion probe does, but the indications are 
there . 

The conclusion to be drawn from these measurements is that the cavi- 
ton is surrounded by steep charge density gradients caused by "digging in” 
of the applied RF field. In our discharge the cavitons are actually 
visible. The structure seen must consist of excited atoms, which form a 
density wall around a dark space in which at least no excited atoms exist. 
The ions and electron are, of course, not visible but the probe measure- 
ment indicates existence of similar density walls as the visible ones. 

The Langmuir probe measurements do not have sufficient spatial resolution 
to prove that both kind of walls are located at the same distance from the 
electrodes. Geometry and pressure of the discharge of reference IV-1 are 
too different for allowing quantitative comparison. 

2. Collimated X-Rays 

The original goal of the effort was to show existence of coherent 
x-rays. We have shown evidence of emission of x-rays in a preferential 
direction. However, this does not necessarily mean these x-rays are 
coherent. 

In the following an attempt is made to explain the mechanism of 
generation of collimated x-rays for the case of the degenerated cavitons. 

Both devices — the high voltage pulsed dc discharge and the low vol- 
tage RF discharge — produced at least collimated x-rays. This in indicated 
by the result shown in Figures III-3 and III-13 through III-16. 

Since Figure I II- 3 is always reproducible, while Figure I I 1-13 through 
III— 15 show a wide spread in phenomena, the discussion on collimated 
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x-rays is limited to the pulsed dc device in the following. 

It is claimed that Figure III-3 could only be obtained by exposure 
to an x-ray source which emits in a preferential direction. 

Figure III-3 was obtained by wrapping an x-ray film around the dis- 
charge tube. It is now first attempted to explain the obtained exposure 
by x-rays generated by impact of electrons on the glass tube. Referring 
to Figure IV-4, assume that the electric field has a configuration which 
forces the free electrons to impact on an area of the glass tube adjacent 
to the negative electrode. Assume further that the field has an infinitely 
steep gradient toward the right side of the impact point "A," which means 
no electrons impacting at the right side of point A and many electrons are 
impacting on the left side of point A. 

In this case, point A can be considered as a point source for isotropic 
x-ray bremsstrahlung. (Only if the impact energies were in the MeV range 
non-isotropic radiation could be expected.) The darkening of the x-ray 
film, at the right side of point A, is of interest. 

Neglecting scattering for the moment and dealing only with absorption 
in the glass tube, the intensity of the bremsstrahlung at point x = 0 
on the x-ray film is 


I = I e 
o 


-kd 


where I q is the intensity emitted at point A, while k is the absorption 
coefficient of the glass and d is the thickness of the glass wall. 

The intensity of the x-ray exposing the film at the right side of 
point A (x>0) is given as a function of x by 


I (x) = I q e 


-k( v4?+d 2 ) 





Fig. IV-4. Bremsstrahlung from Impact on Glass Wall. 
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The sensitivity of the photographic material (in the "straight" portion 
of the sensitivity curve) is logarithmically related to the energy expos- 


ing the film. 

Therefore, the ratio between number of exposed grains (density) on 
the left side to the number of those on the right side is at best 

log _ _ ky'(x 2 +d 2 ) 

o 

which is a linear dependence and not at all the steep density gradient 
which is observed. Scattering in the glass was neglected, as well as 
the fact that the points left of point A are also isotropic x-ray sources 
and will add to the exposure of the x-ray film at the right side of point 
A. This will even further reduce the expected density gradient on the 
film. 

The conclusion is that the exposure shown in Figure III-3 cannot 
possibly be caused by bremsstrahlung generated by impact of electrons on 
the glass wall. 

However, the figure can be explained at least qualitatively by the 
assumption that a "degenerated caviton" existed under the conditions at 
which the exposure was obtained. Figure IV-5 shows a number of such de- 
generated cavitons for the RF discharge, where they are visible. This 
figure has to be compared to Figure IV-6 (a repetition of III-4) which is an 
x-ray exposure- By comparing these two figures, it is readily seen that 
such degenerated cavitons must be responsible for the observed collimated 
x-ray radiation. 

A possible mechanism could be bremsstrahlung (or cyclotron radiation) 
which is caused by electrons vibrating back and forth across the boundary 
line of the caviton. In order to produce x-rays under such conditions. 
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Figrue IV-6 Three dimensional sketch of X-ray pattern around the tube 
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we would have to assume very strong field gradients at these boundaries. 
But, if such strong gradients exist, the generated x-rays would have to be 
collimated and even coherent to some extent, since all participating elec- 
trons are vibrating parallel to each other. 



V. DESCRIPTION OF THE STATIONARY DISK-SHAPED ENRAD 
RESONANCE-PATTERNS IN TERMS OF CAVITONS 


1 . Derivation of the Fundamental Nonlinear Equation 


In order to analyze the developed caviton, its stationary states, and 
their stability, we start from a set of three equations: 

36n 


-sr 5 - + V [ (n + fin )v ] - 0 
dt o o e 


( 1 . 1 ) 


3v T 

(n + 6n ) r-“ + (n + 6n ) + VSn = 0 

o o dt o o m m e 


( 1 . 2 ) 


V* E = - 47re 6n (1-3) 

e 

Here we have divided the electron concentration m into three components 

n = n + 6n + 6n , (1.4) 

o o e 

where n Q is the average concentration of electrons balanced by the average 
ionic concentration in the discharge tube, 6n Q is a very s low fluctuation, 
neutralized by an identical fluctuation in the ionic concentration (also in- 
cluding a large time- independent part defining the caviton) , and 6n^ is the 
fast fluctuation which cannot be followed by the ions. We have considered 
both 6n^ and the velocity v^ of the electrons small and we have neglected 
terms quadratic in them. The first equation is the equation of continuity 
for the fast (jitter) motion, the second is the equation of motion expressing 
the momentum balance, with a field (drift) term and a pressure (diffusion) 
term included in addition to the linearized inertial term (T^ in energy 
units). Equation 1.3 is Poisson's equation. We neglect 3Sn/3t, i.e., 
consider 6n time- independent. 

We take the divergence of equation 1.2 and subtract the time -derivative 
of equation 1.1: 
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9 2 6n T 

- - + - V* [ (n + fin)E] +^V Sn = 0 

2 mo me 

o t 


(1-5) 


Using equation 1.3 we obtain: 


92^ . 2 T 

— t + in®_ (n + 5n) E - — V(V-E) 
3t 2 m 


= 0 


( 1 . 6 ) 


A solenoidal part of the rectangular bracket does not make sense physically 
in our case. Neither does a harmonic potential part. Therefore we are 
looking for a solution satisfying the equation: 




(n + 6n) E 
o 


T 

— V(V-E) 
m 


0 


(1.7) 


We separate out the fast part of the variation in E(t): 


E = 77 E (t) e i “ pt + C.C., (1.8) 

Z o 

where E (t) is a slowly varying amplitude with a negligible second deriva- 
o 

tive implied, and c.c. is the complex conjugate of the preceding term. The 
average ponderomotive force of the electric field ^(t) on the plasma of elec- 
trons and ions can be described by an average potential energy'** (see also 
section 7 ) . 


U = e E 


A 2 
4mw 


2m w 


<e 2 > , 


(1-9) 


where the average <> is a time average. Denoting by and the ionic 
and electronic temperatures in energy units, we expect a Boltzmann distri- 
bution of the plasma density 

, ~ -U/(Ti + Te) ,, _ 

n + on = n e (1.10) 

o o 


It is convenient to introduce the dimensionless field amplitude 

A = eE* / [2 wvfr(T " + T ) ] . We obtain 
o l e 
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2i (a) 


3a 

3t 


2 + 2 - I A I 2 t 

0) A + W e 1 1 A 

P P 


£ 2 u ) 2 V(V*a) 

D p 


0 . 


( 1 . 11 ) 


2 1/2 

The electron plasma frequency U)^ = (47ine /m) and the Debye 


= /T e / (47fne ) define the dimensionless variables T = to^t/2 
This yields the fundamental equation 


length 

and p = r /^ D - 


3 A - 1 A I 2 ^ 

i - 5 “ + V p (V p -A) + (1 - e ) A = 0 (1.12) 

the stationary (steady state) solutions of this equation will be analyzed 
in the next section. 

This equation is known as the nonlinear Schrodinger equation, which 
describes both the temporal evaluation and the stationary states of normal- 
ized field amplitude A. 

As we shall see in the next section. Equation 1.12 admits stationary 
localized solutions similar to bound stated encountered in quantum mechanics. 
We will call these stationary localized states of the electric field "cavi- 
tons. " 


2 . Stationary Solutions of the F un damental Equation: Cavitons as Dynamic 

Equilibrium Configurations 


Consider first the one-dimensional case 


. 3a 
1 3 1 


+ 


3 2 A -I a| 

3"T + (1 " e 


) A = 0. 


( 2 . 1 ) 


z 

This equation admits stationary solutions in which the dynamic pressure of 
the trapped electromagnetic field is balanced by a variation of the electron 
thermal pressure. These solutions, known as cavitons, are a special case of 
solitons, and have the form 


A(z,T) 


V(?)e 


ikz + icoT 


£ = z - 2kT. 


( 2 . 2 ) 
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Equation 2.1 becomes now 

2 2 

+ (1 - e~ - w) y = 0; w = 0 ) + k 2 . (2.3) 

ar 

2 

For the case of the small -amplitude (u «1) Langmuir soliton (caviton), we 
expand the exponential in series. We obtain in lowest order an equation 
studied by Zaharov** - , which presents a cubic nonlinearity 

h 2 

^=r + y « wy. (2.4) 

as 2 

This equation admits an exact solution 

y(S) = C sech (|c|5//2); w = |c| 2 /2, (2.5) 

where C is an arbitrary complex constant. Indeed, equation 2.4 is satis- 
fied if we substitute this "soliton" solution. Note that for any value of 
w in equation 2.3 or 2.4 only one solution is possible with the form seen 
in Figure 2.1. As we shall see, the situation is different, and more 
interesting in three dimensions, where the spherically symmetric solutions 
which are regular (finite) in the origin constitute a discrete set, with 
quantized radial slopes of the allowed solutions at r = o. Spatial depend- 
ence of the normalized field amplitude of the one-dimensional caviton. The 
dotted line is the low amplitude limit C sech (C?//2. The full line is 
obtained by numerical integration with the use of the Runge-Kutta proce- 
dure. Spatial variation of the plasma density for |c| 2 = 2 in Figure 2.1. 

The spatial dependence of the plasma density in the one-dimensional cavi- 
2 

ton is calculated in stationary conditions by using equation 1.10. The 
results are plotted in Figure 2.2 for the large-amplitude caviton with 

i 1 2 

| C | =2 whose field is represented in Figure 2.1. Note that only 13.5% 

of the plasma density is left in the center of the caviton. For higher 
amplitude solutions this percentage will further decrease exponentially. 
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according the equation 1.10 and 1.9, until a good vacuum is left in the 
caviton. 

We consider now the three-dimensional, spherically symmetric Langmuir 
caviton. Equation 1.12 becomes 


3a 1 3 2 3a 

3t 2 3p 3 * 
P P 


2 - 1 a| 2 

— A + (1-e |AI ) A = 0 

P 2 


With the boundary conditions : A bounded for p = 0 and A 0 for P °°, one 

obtains a discrete set of solutions with n = 1,2,3... modes or zeros of 
A between p = 0 and p =°°. 

The result of a numerical calculation 3 for the fundamental mode and 
the first "excited" mode are represented in Figure 2.3. The electric field 
points in the radial direction but vanishes in the origin and, for the case 
of nonfundamental caviton states, also on the nodal spheres. 

For prx» the solution in Figure 2.3 decreases like spherical Bessel 
functions of order unity. The characteristic length of the high- field 
region is 


L = £ 0) /u) = ( 

D p 


4TTn e 
o 


47Tn e 1/2 

— f— > of 1 




This means that the size of the caviton is independent of the surrounding 

~ —16 —13 

plasma density n^. For the electrons T = 2000°K*1 . 38 . 10 erg = 2.76 10 

-27 7 -1 

erg = 0.1725eV and with m = 0.9.10 g we obtain 1.75.10 cm/s w . This 
yields L = zcm for a frequency of 1.34 MHz which corresponds to the largest 
neutron output of ENRAD. The observed width of the dark space in ENRAD 
is approximately 1 cm, i.e., of the same order of magnitude as a Langmuir 
caviton. 


The antinodal regions in Figure 2.3 correspond to spherical plasma- 
depletion layers similar to Figure 2.2, and the nodal spheres contain the 
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normal plasma concentration n Q . The similarity of spherical Langmuir soli- 

ton states to electronic states in atoms is particularly interesting. 

1/2 

The normalized electric field (3/2) E vs. p= 03r/w^ D , of the ground state 

and first excited electrostatic mode with spherical symmetry, and with 

a ) 2 /(jl) 2 = 1/2. The electric field is radial; 3 2 = e 2 /k7Tm 2 (T + T. ) . 
p e e l 
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3 . Electromagnetic Field Configuration and Plasma Generation in ENRAD 


At the ENRAD operating frequencies of 0.5-1. 7 MHz, the vacuum elec- 
tromagnetic wavelength is 6 00m- 17 6m. The distance between the ENRAD 
electrodes is between 5 and 50 cm, which is much smaller. Therefore, we 
can approximate the electric field configuration between the electrodes 
with a statis electric field. 

Let the axis of the discharge tube be the z axis and the axial point 
0, equidistant to the electrodes, be the origin of a system of rectangular 
(x,y,x) and cylindrical (r,6,z) coordinates as shown in Figure 3.1. 

First we obtain the electric field E in the vicinity of the origin 0 
in the median plane z = 0 by summing the contributions of corresponding 
elements of arc dl* and dl" on the cir cumber ence of the two electrodes of 
radius a = 10 cm, carrying a charge O per unit length. 

dl 


e = a 


dl 


-v -> 
r-r 1 


r-r" 


i 3 

r-r" 


r-r" 


= - OR 


2 ->-2 | v 3/2 

r' +r -2r.r ' ) 


E = 


aaR 


<»W/4> V2 


^ 2 2 2 2 * 
1 , 3 r -2racos0 . „ _ a r cos 0 

flu I -L _ _ T / • J 

2 2.2,. . 2.2 ,.,2 


a +R /4 


(a+R/4)J 


E = 


- 2iTaaR 


la 2 t R 2 /4> 3/2 


1 t 


pr 


2 2 
a +R /4 


„ - 3.75a 

P = 


2 2 
a +R /4 


1 3 I j 3 -*■ 

Here we have used the equality |r-r' | - |r-r" | valid for r in the median 

plane, and r'-r" = R. The coefficient o can be calculated for the distances 
5 cm < R < 50 cm used in ENRAD. We obtain 


R 

5 cm 

a ~ 24.5 cm 

50 cm 

P 

2.03 

0 

-0.983 


Consequently, the field distribution in the median plane z = 0 is charac- 
terized by a maximum in the origin for R > a^6 = 24.5 cm and by the pres- 
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ence of a maximal circle with a minimum in the center 0 for R > a/6. 

For R = a/6 the field maximum in 0 is flat and of higher order, since all 
second derivatives of the field vanish. 

The dependence of E(r,z) on the cylindrical coordinate z and the 
radial field component can be determined in the vicinity of the origin 
by performing a more complete expansion of the first form in equation 3.1 
or by using VxE = 0 and V*E = 0 with equation 3.3. Both methods yield, 
to second order. 


E = - 
z 


27Ta<JR 


, 2 2 ... 3/2 

(a +R /4) 


1 +p 


2 o 2 

r -2z 

2 2 
a +R /4 


47TaQR 

E r " ,2 2.5/2 


prz 


(a +R /4) ' 

In three dimensions, the maximum observed in 0 for R > a/6 turns out to 

->2 

be a saddle point, both for and for E . The same holds true for the 

minimum observed in 0 on the z = 0 plane. 

The main features of the whole field configuration are represented 

qualitatively in axial cross section on Figure 3.2 for the three cases 

considered. The maximal-field points obtained in any plane z = const 

generate a surface in three dimensional space. The cross section of this 

surface is represented in Figure 3.2 by full lines, and the median 

|e| = eI = const, surfaces are dotted. The arrows indicate the di- 

1 1 l z=0,max 

rection in which |e| increases on these lines which are the geometrical 
locus of the points of maximal electric field lines. For R a/6 the locus 
of the maximal- field points is a cylindrically symmetric surface which 
resembles a rotational hyperboloid with one sheet. For R = a/6 the neck 
width of the rotational surface becomes zero and for R a/6^ the locus consists 
of two distinct funnel- shaped surfaces connected by a line on their axis. 
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From the shape of the |e| = max and |e| = const surfaces, we note 

that in the center of the ENRAD device there is an extended, quasi-homo- 

geneous high- field region, separated from wall and electrode surfaces 

which are well known plasma recombination regions- Therefore, after the 

application of the high-frequency field in ENRAD, the highest plasma 

density is expected in the median (z = 0) plane. 

2 1/2 

The plasma frequency 0) = (4TTne /m) increases together with the 

electron concentration n in every point after the high frequency field 

has been applied and after the first breakdown has occurred. In the 

process of plasma build-up, the increasing plasma frequency values 

U) (n) will first approach the value to of the applied field in some point 
p o 

P on the median plane z = 0. As soon as to^ approaches to^ the dielectric 
constant of the plasma vanishes and since EE = const, the electric field 
in the point P is enhanced (mathematically it approaches infinity) . 

At this moment, a resonance between the electromagnetic wave and the 
electrostatic plasma oscillations characterized by to^ occurs in the vicin- 
ity of the point P. As we shall see in the next section, this resonance 
leads to a dramatic enhancement of the electric field in the vicinity of 
P and to the process of linear conversion of electromagnetic waves into 
electrostatic (Langmuir) waves. 

4 . Electrostatic Resonance, Linear Conversion from Electromagnetic to 

Electrostatic Waves in ENRAD, and the Development of Localized Fields 

In order to better understand the resonance process at the point P 
defined in the previous section, in which to^ first approaches to, we first 
review the elementary derivation of the electromagnetic wave dispersion 
relation and plasma dielectric constant. By applying the curl and a time 
derivative, respectively, to the two Maxwell equations 






%. • 
V 
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V X E = - — 


c 3 t ' 


„ £ 1 3E ^ 4tt -j- 

VxB = caE + r ] 


we eliminate the Vx 3B/3t term and obtain 


7 x (V x ^ - ~ I?- 


2 3 t 


2 3 t 


Taking 3 = nev and using the equation of motion 


(4.1) 

(4.2) 


dv e 

dt “ m 


we obtain with VE = 0 

_2-* 1 3 2 e 4 Tine 2 -*• 

* 2 » 2 2 E 
c dt me 

Substituting E = E^e i (k r - OJt) we obtain the dispersion relation 

2 2 2 2 
0 ) = 0 ) + c k 

P 


and the phase velocity 

2 

V ph 


0) 


2 / 2 
1-0) /(A) 

P 


_ 

e 


, 1/2 


This identifies the index of refraction e and the dielectric 
constant of the plasma 


2 2 

e = 1 - o) /o) . 

P 

2 

Similar results obtained by replacing c with 3kT/m in equation 4.5 
and 4.7, holds for Langmuir waves (see Section 6 ). Consider the ENRAD 
discharge right after the high frequency field has been applied, when the 
density n of electrons (and ions) starts to increase and the plasma fre- 
quency in the z = 0 plane approaches the frequency (a) q of the applied field. 

->• -v -v 

Along the discharge tube we have roughly E = 




The electric field configuration is close to a static configuration with 
the field lines directed along the axis, between the electrodes. Assuming 
that a)p approaches tO o in a cross-sectional plane close to the median, we 
expect an enhancement of the electric field 


E 


D 

e 


p 

2 2 
i - a) /or 
p 


2 2 

in this layer. The enhancement factor is 1 - 0) /a) . This enhanced field 

P 

will further increase the degree of ionization and consequently 0)^ will 
reach the value ai Q . Note that the energy density ED/8 tt is also increased 
in the resonance layer. Since the amplitude of D is constant, i.e., the 
same as outside the layer, the swelling in the energy density is propor- 
tional to that observed in the electric field amplitude. 

The increased electric field in the resonant layer will excite elec- 

—>■ 

trostatic (Langmuir) waves. Indeed, the enhanced field E will accelerate 

—y 

electrons in the -E direction, creating a compression region (accumulation) 
on the -E side of the resonant layer and a depletion region on the +E side. 
Due to the electrostatic repulsion between the compressed electrons on the 

. -4. 

-E side, and due to the pull of the uncompensated positive ions on the +E 

side, the electrons will rebound after T /2 = TT/u) , i.e., after a half- 

P P 

period of the Langmuir wave. Indeed, this electrostatic push-pull is the 
restoring force of the electrostatic waves in a plasma. Due to the equality 
w = the applied electric field will have changed its direction from 
E to -E in the same time. Therefore, the electrons will be accelerated 
again during their rebound. The enhanced electric field thus continuously 
builds up the plasma oscillations by acting on the electrons with the right 
phase at resonance. The electrostatic oscillation would increase 
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indefinitely in amplitude, like a lossless resonator, if there were 
no energy-loss processes and non-linear phenomena saturating the ampli- 
tude at some very high values, e.g., 2 or 3, orders of magnitude above 
the applied electric field amplitude level. Such processes are: 

A. Generation of Langmuir waves which propagate away from the 
resonant layer in both directions. The coupling of the electromagnetic 
field to these waves is given by the excitation of electrostatic waves 
by the enhanced electric field which is in the resonant layer, part both 
of the electromagnetic and electrostatic wave. This excitation, described 
in detail above is limited to the resonant layer, because the resonance 

condition a) - a) is only satisfied in this layer. The coupling is linear, 
p o 

and therefore we speak about linear conversion from EM to ES waves. The 

Electrostatic waves generated in this layer propagate down the plasma 

density gradient adjoining the layer in both axial directions, with in- 

2 2 ly ^ 2 

creasing group velocity v ges = ^ se I 1 - W^(z)/o> ] , decreasing energy 

density by the same factor, and decreasing wavelength A = (2^/c ) 

- 1/2 Se 
2 2 ' - 1/2 

[1 - W (z)/W ] , where c - (3kT /m ) is the effective electronic 

p se e e 

sound velocity. The electromagnetic energy streams into the resonant 
layer along the field lines of the Poynting vector (c/4TT)ExB at the speed 
of light and leaves the layer at the much lower speed v^ of the order of 
the electronic thermal velocity, in form of Langmuir waves. Energy flux 
conservation requires then that the electrostatic energy density in the 
Langmuir wave must be larger than the electromagnetic energy density by 
the ratio of these speeds, i.e., 

V C (1-W 2 /0J 2 ) 1 / 2 

' gem — = c/c ~ 3 . 10 10 cm/s/3 • 10 7 cm/s = 10 3 

V ges c a- w / u> ) / 
se p 

This energy enhancement is equivalent to an additional increase of the 
electric field amplitude by 1~- orders of magnitude, actually only a 
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2 2 -I 

factor of 10 in addition to the enhancement due to the (1 -oj /a) ) factor. 

P 

The physical reason for this increase in the energy density is the kinetic 
energy of the electrons in the Langmuir wave and the electrostatic energy 
of the charge density accumulations in the Langmuir wave, both contribut- 
ing equal amounts in average. The actual factor is reduced to 10 due to 
the competition of other energy loss process, listed below. 

B. Charged particle acceleration - to be discussed in detail in 

4 

Section 8, leads to transit damping : suprathermal electrons are accele- 

rated by passing through the resonant layer with the right phase. In 
addition, as we shall see in the next section, the ponderomotive force of 
the enhanced electric field acts on the plasma and digs a cavity leading 
to the creation of an electrostatic caviton. This large density pertur- 
bation corresponds to a localized ponderomotive potential (see section 7) 
which also contributes to the acceleration of electrons and ions observed 
experimentally in ENRAD. Fast electrons and ions carry out about equal 
amount of power, as we shall see. 

C. Resistive (ohmic) heating in the walls of the density cavity. 

D. Generation of ion-acoustic waves and of other oscillations of the 

plasma. These processes are based on the coupling between various resonant 

modes of the plasma. If two modes are involved, of frequencies and 

and wave numbers k. , and k_ , the relations a) = at + (jlL and k = k, + k 
12 o 1 2 12 

must be satisfied. 

The enhanced electric field present in the resonant layer initiates 

4 

the linear conversion from EM to ES waves on a short time scale of the 
order of 10 periods of the electronic plasma frequency to^, determined 
essentially by the inertia of the electrons. The subsequent development 
of a density minimum in the resonant layer is about two orders of magni- 
tude slower, i.e., of the order of 10 periods of the ion plasma frequency 
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2 1/2 

to . = (4 n.e /m.) , where n is the ion concentration and m their mass. 

P3 i i 

5. Formation of Cavitons in the Nonlinear ENRAD Regime 

The strongly enhanced electric field present in the resonant region 

~ 2 

with to - to exerts a ponderomotive force f “ - V <e >/87T per unit volume 
P o 

on the plasma, thereby creating a density cavity in the resonant layer. 

This density cavity is favorable to the further trapping of electromagnetic 

field, reducing the convective loss of electromagnetic and Langmuir wave 

energy. This trapping results in an even higher field inside the cavity, 

which in turn creates a deeper cavity. This process continues in a self- 

consistent way until saturation caused by nonlinear phenomena sets in. 

This final stationary state obtained this way is a combined nonlinear 

field and plasma state known as caviton, a particular case of solitons 

well known in many nonlinear physics problems. As we have seen in section 

2, the electric field inside the caviton is very similar to the quantum- 

mechanical wave function defined inside a potential well which corresponds 

to the density cavity in our case. The splitting of a single caviton 

with formation of two cavitons when the pressure is slowly raised in the 

ENRAD discharge is similar to the transition from the ground state to the 

first excited state of a sequence of energy levels in quantum mechanics. 

As we have seen on the example of the numerically calculated spherical 

caviton, this feature is obtained also theoretically on the basis of the 

nonlinear "Schrodinger" equation derived in section 1. The nonlinearity 

th 

is due to the interaction of the field with the plasma. The n excited 
state has n nodal spheres. 

In the case of the ENRAD device, the ground state corresponds to one 
dark space in the z = 0 plane bounded by a pair of luminous disks. The 



104 


dark space is the high field region in the caviton and the two disks are 
the walls of the density cavity, made luminous by the excitation of atoms 
under the influence of mainly electronic bombardment and also ionic bom- 
bardment caused by the high electric field. In fact the luminous disks 
arise from the plasma distribution partially penetrating into the high 
field region and adjacent plasma sheets extending into the background 
gas as in Figure 5-1. In the high field region, inside the caviton, the 
concentration of exicted ions and neutral atoms is very low, and there- 
fore this region appears dark to the eye. 

The excited state of number n represents a sequence of n high-field 
regions in ENRAD, each of them limited by two luminous disks. 

Figure 5.1 is disk-shaped plasma cavitons obtained in the ENRAD 
device: a) distribution of the squared electric field amplitude; b) 

distribution of the total plasma and gas pressure, p = n^kT^ + n^kT_^ + 

n D kT D + n D +kT D + N D kT D ? distribution of the light output in the 

blue and red spectrum areas, and the total luminosity; and d) distribu- 

tion of the total plasma concentration n = 1/2 (n +n.+n +) . 

e l D ' 










106 


" 


6. Derivation of the Electrostatic (Langmuir) Waves 

We have repeatedly used the concept of electrostatic waves and the 
expressions of their group velocity and phase velocity in the present 
chapter- The plasma caviton comprises, in essence, a standing, localized 
wave of this kind. Therefore, we include here a brief derivation of elec- 
trostatic waves in plasmas, which may be skipped by the reader who is 
familiar with the subject. 

Starting with the generalized Ohm's law 


_5L|1 = | + v x | + l_ Vp ___l_ (e.D 


ne 


2 8t 


Vp 

en e cne 


where j is the current density and T) is the resistivity, we take B = 
0, and consider small, fast, fluctuations of the electron pressure 

2 . 


<5p e = Y kT e VSn = - 3 ^ Vp = - J ‘Vp, 


( 6 . 2 ) 


2 _ 

where v^,_ = 3 k Te/m, and the adiabatic coefficient y = (s/z + l)/(s/2) = 3. 

th 

The number of degrees of freedom s has to be taken = 1 since the fluctuations 
are faster than the collision time and the motion is, therefore, adiabatic 
and one-dimensional. The fast fluctuation p of the electronic charge 
density satisfied the continuity equation for the fast fluctuations 


9£ 

at 


+ 


V(6j) = 0 


(6.3) 


allows us to eliminate 6j from the variation of Equation (6.1) by applying 
the divergence: 


a 2 p 2 2„2 

37 + W p P ' V th V P+V 


f£> = 0 
at 


(6.4) 


2 2 2 

Here V is the frequency of collisions V = ne /m, and = 47fne /m is the 
squared plasma frequency. Note that Equation (6.4) differs from Equation 
(1.5) only by the additional presence of the dissipative collisional term 
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with V. The corresponding dispersion relation for plane wave solutions 
substituted into Equation (6.4) is 

(d 2 = 0) 2 + v, , 2 k 2 - iidV (6.5) 

p th 


Therefore, if we neglect the small dissipative (damping) term are obtain 
the phase velocity 

2 _ “p 2 . V_ V th 

V ph (e } .2 2 V th " 1V ph k 2 2 “ 'th 

k ;k 1-0) A) 


2 . 
v 4.u /£/ 


( 6 . 6 ) 


where £ is the dielectric constant defined in Equation (4.8) . The group 


velocity v is given by 
v 2 k 

v (e ) — 

g 3k o) 2 2 

p + V th k 


2 2 

k 2 V th , “p 

m V th ° “7 = v th 1 - ~ = v th^ 

v , (e ) 0) 


ph 


(6.7) 


Note that Equations (6.5)- (6-7) turn into the corresponding equations 

(4.6) and (4.7) applicable for electromagnetic waves, if we replace v , 

th 

by c, the speed of light in vacuum. In particular, the group velocity of 
electomagnetic waves is 

2 U) 2 

v (em) = c — = — — — = c 1 - ~^r - c /e (6.8) 

5 “ v ph (em) U) 2 


7. Derivation of the Ponderomotive Potential 

The concept of the ponderomotive potential used in the previous sec- 
tions arises in a natural way, if we consider separately the fast motion 
of the electrons at frequencies of the order of the plasma frequency, and 
the slow motion in which both electrons and ions participate. 

The high frequency (jitter) velocity of the electrons is given by 
the equation of motion 




e 


7 # 

■&' 


(e/m)E 


(7.1) 
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which produces an additional meansquare velocity 


<v 2 > = (e/mw) 2 <E 2 >. 

e 


(7.2) 


This corresponds to an additional pressure contribution p ff 

p = (1/2) n m<v 2 > = (1/2) (n e 2 /m03 2 ) <E 2 > = 

eff e e e ott 


(7.3) 


In turn, this effective pressure causes a force -V P e ^ f P er unit volume. 
If we divide by the concentration of electrons n, we obtain the average 
force per electron 


^ V R <E^ > <E^ > 

2mOJ e 


(7.4) 


where v is the slow component of the electron (and ion) velocity in which 
o 

the jitter has been averaged out. The form with n in the denominator is 

2 2 2 

applicable in the case of a plasma caviton, for which 0) = (0^ = ^Trn^e /m. 

The potential energy U from which the ponderomotive force can be derived 
is thus identical to equation (1.9) 


_ e 2 <E 2 > 

u = — r 

2mW 


<e 2 > 

87Tn 


(7.5) 


also known as the slow pseudopotential, felt only by the electrons. The 

U/kT e 


electrons tend to assume a Boltzmann distribution n = n e 

e o 


e , but the 

ions pull them back partially, through the amibpolar potential <J) and a 
corresponding ambipolar electric field which always appears if elec- 

trons and ions are acted upon with different forces. The ambipolar 
potential <j> restores electric quasineutrality, and we obtain a distribu- 
tion of the electron and ion concentrations 

U-e<|> ~ 

- - ■ _ m ■ = n. = n e kT. 
n=nekT i o i 
e o e 


(7.6) 



109 


From the equality of the two exponents in Equation (7.6) we infer that 

et = St , _S_ . (7.7, 

T T. T +T. 

e l e l 

This yields the final form of the spatial dependence of the carrier con- 
centration 


n = n . - n e k (T +T . ) 
e l o e i 


(7.8) 


which has been used in Equation (1.10). 




VI. SUMMARY AND CONCLUSIONS 


During the effort reported herein, investigation concerning the 
properties of a radiofrequency discharge of a special configuration, 
called ENRAD, were made. The results obtained were compared to measure- 
ments previously obtained with a pulsed high voltage discharge employing 
the same conf iguration. 

The observations made during these investigations are summarized as 
follows: 

1. Emission of neutrons by the RF discharge was observed. 

2. Indications of some weak x-ray emission were also obtained. 

3. A visible plasma structure consisting of pairs of disks separated 
by a dark sheath was observed. 

A model to explain the observed phenomena was proposed. This model de- 
scribes formation of plasma cavitons caused by "digging in" of the applied 
RF field. The condition for this to occur is that the applied frequency 
is equal to the plasma frequency. 

The caviton is a plasma region of drastically reduced density sur- 
rounded by steep gradients in electric field and charge carrier density. 
Charged particles trapped inside the caviton can obtain high energies due 
to these field gradients. They may escape from the caviton and bombard 
the surrounding high particle density walls. The resulting collisions 
are sufficiently energetic to cause nuclear reactions to take place. 

It is concluded that one of these reactions is the DD reaction causing 
neutrons to be emitted. 

There is some experimental evidence these nuclear reactions take 
place, especially during the collapse of the cavitons. This is con- 
cluded from the fact that introduction of small concentrations of nitro- 
gen into the discharge increase the neutron yield. 



It is also concluded that in principle these cavitons are capable 
of emitting collimated, maybe even coherent, x-rays. This could be 
especially shown for cavitons formed at higher voltages. 

The conclusion is that the ENRAD configuration is a plasma contain 
ment which is capable of producing nuclear reactions at high particle 
densities and low temperatures. The applied electric fields are modest 
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A series of radioactivity measurements have been performed for a 
number of high purity gold and indium foils which have been placed into 
positions around the deuterium gas discharge device for varying periods 
of time while the device was operated. The objective of the measurements 
was to monitor the device for neutron emission, which had been reported 
by others; having been deduced from evidence attributed to neutron inter- 
actions in detector media or radioactivation of metallic detector foils or 
disks (coins) which were subsequently counted. 

The detector system used in the measurements, as shown in block diagram 
form in Figure 1, consisted of well-shielded, high efficiency beta-gamma 
scintillation detectors with attached instrumentation combined into a beta- 
gamma absolute coincidence spectrometer. The beta detector was CaF^ scin- 
tillation crystal, while the gamma detector used was a Nal(Tl) scintillation 
crystal. Beta and gamma ray singles counts (Ng, N^) and beta-gamma coinci- 
dence counts (or counting rates) were measured with single channel and 

coincidence analyzers, while at the same time, the ungated and gated gamma 
ray pulse height (energy) distributors were measured with a multichannel 
pulse height analyzer. 

The beta-gamma coincidence system as described, and shown in the figure, 

198 

was selected to provide a high degree of selectivity for Au, through its 

100% branching through the beta-gamma cascade and the systems ability to 

198 

probide a direct value for the Au source strength N^, i.e., Ng» N^/N^ = 

n d. 

Radioactivation of pure indium and gold, such as the foils used, would 
form two general types of radionuclides, depending on the energy group 
characterizing the neutrons; one for thermal and the other for fast neutrons. 
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The thermal neutron group, in general, form y or negatron decay species, 
through the (n, y) reaction, while the fast neutron group form either 
electron capture (EC) or positron decay radionuclides by the (n, 2n) reaction. 

The principal (n, y) species formed for In and Au are the 54-minute 
116m 198 

in and the 2.7-day Au radionuclides. The former emits several gamma 
rays following its negatron decay with energies (% branching) ap 4.7 (36%), 

819 (17%)m 1090 (53%), 1293 (80%) keV, while the latter exhibits only one 


gamma ray at 411 keV (100%) . 

The (n, 2n) fast neutron induced reaction products for In would be 
derived from 

(4.23%) 113 In(n, 2n) 112m in(20. 7 min) IT(156 keVy) 

112 In(14.4 min., 6~(44%) , 8 + (22%) , EC (34%) 

511 keVy (44%) 

617 keVy (6%) 


and 

(95.77%) 115 In(n, 2n) 114m ln (50-day) IT (96. 5%) (191.6 keVy's) 

EC (3. 5%) 724, (58 keVy's) 

114 In(72 sec.) g" (98%), B + (.004%), 1299 keV (.2%) 
Gold (n, 2n) reaction products are those of ^^Au and ^^ m Au: 

(100%) 197 Au (n, 2n) 196m Au (9.7hr) IT (175 keVy) 


196 Au (6.1 day) 8" (6%) -426 keV y-ray EC (93%) 688.7 
(26%), 355.7 (67%) keVy's 

8 + (5 x 10 _5 %) 


A common component introduced into gold coins is copper, which forms 
64 4- 

and (n, y) 12.8 hour Cu positron activity (3 19%) by means of thermal 

neutron absorption. The fast neutron group would also produce the same 

62 

positron activity, but would, in addition produce 9.76-minute Cu, i.e., 

(30. 9%) 65 Cu (n, 2n) 64 Cu (12.8 hr.) g“(38%) - no gamma rays; EC (43%) - 

1340 keVy (.5%); B + (19%) - 511 keV 
annihilation gamma. 
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(69.1%) 63 Cu (n, 2n) 62 Cu(9.76 min.) g + (97%) - 511 keV annihilation 

gamma . 

Thus, in general, positron annihilation observed in the decay of pure 
indium is possible but for gold it would have to be due to activated im- 
purities, one of which could be copper, which could be produced with either 

62 

thermal or fast neutrons. The shorter lived copper radionuclide Cu(9.7m) 
would be indicative of fast neutron activation since it can be produced only 
by the fast neutron (n, 2n) reaction. 

The first radioactivity measurements were made with a gold coin, 

several gold foils and several indium neutron detector foils, which had 

been placed around and near the discharge device for approximately 12 hours 

of operation. The gold, having a lifetime of 2.7 days, would allow for 

longer counting times , while the indium foils with a lifetime of 54 minutes, 

should provide higher counting rates for shorter activation periods. 

For the first measurements the discriminator for the CaF^ beta ray 

198 

detector was adjusted high enough to exclude output pulses due to Au 

gamma ray interations in the thin CaF^ detector. This assured only beta- 

gamma gold coincidence counts could occur. Even though a low level back- 
137 

ground due to Cs gamma rays was observed in the ungated singles gamma 

spectrum, the gated (coincidence) spectrum, however, was essentially free 

137 

of effects due to the Cs gamma rays. 

A series of background counts for comparison were obtained prior, 

during, and after the initial measuremtns. The Au and In foils were counted 

before and after their exposure to the operating deuterium discharge. After 

a delay due to cross-town transport, the In foils were counted first due to 

198 

their shorter lifetime (54 min.), with the longer lived Au (2.7 days) 
being held for later counting. The first counts were for 10-minutes, how- 
ever the total counts obtained were marginal for satisfactory statistical 
evaluation. To improve counting statistics, the In foils were combined (as 



were the Au foils later) and the counting interval was extended to 40 minutes 
The counting schedule is summarized in Figure 2. The obtained results are 
plotted in Figures 3, 4, and 5, along with the background counting results. 

The first of the indium foils counted. In #002, exhibited counting 
rates approximately equal to the average values of the background. However, 
the next two foils counted, numbers 001 and 003, both exhibited counting 
rates considerably (significantly) above the background average values. 

Also, when these foils were combined with #002 and counted for a period of 
40 minutes, the counting rates observed were again above the considerably 
average background values. 

When the smoothed (3-point) spectra, obtained simultaneously with the 

counting rate data, as well as that of the background are viewed, each is 

observed to exhibit spectral features at approximately channels 38 and 75, 

where the latter peaked feature has previously been shown to be due to the 

137 

662 keV gamma ray from a Cs contribution of the counter-shield system. 

137 

(See Figure 6) Based on using the Cs peak as an internal standard, the 

centroid, or central channel of the peak distribution, of the other peaked 

features is shown to fall at approximately 340 keV. In addition to these 

features, which appeared in the background spectrum, indium foils numbered 

001 and 003 exhibited additional features centered at approximately 

channels 46-48 and 55-56, respectively, while the combined foil 40-minute 

spectrum exhibited the peak at channels 46-48. (See Figure 8) 

A separate measurement performed with an indium foil which had been 

activated by exposure to a thermalized 5 Ci Pu-Be source yielded a peak 

centroid for the 417 keV gamma ray of ^^^ m ln(54 min.) at channel 46. A 

22 

gamma ray spectrum acquired with a Na source yielded a peak centroid for 
the 511 keV positron annihilation peak at approximately channels 74-76. 

(See Figure 8) 
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As summarized in Table I, it was observed that indium foil # 002 
exhibited neither counting rates greatly different from background nor 
spectral features not found present in the background gamma ray spectrum. 
However, indium foils numbered 001 and 003, whose counting rates were 
significatly greater than background also exhibited spectral features that 
were not found present in the background. Foil # 001 exhibited two such 
features, one in the channel range of the 417 gamma ray peak of ^^^Xn and 
the other in the range of the 511 keV positron annihilation peak. Yet, 
only the latter appreared in the spectrum of foil // 003, while only the 
former appeared distinctly in the spectrum for the combined foils. Further, 
it should also be noted that the 54 minute ^^ m in activity is produced 
predominantly by thermal and/or epithermal energy neutrons, while positron 
activity is most frequently associated with charged particle or fast neutron 
((n, 2n)) induced reactions or high energy gamma ray induced electron- 
positron pair formation. 

The counting rate obtained for the 40-minute measurement of the gold 

coin (which was started after two (2) hours and eighteen (18) minutes waiting 

time) was 2.4, 1.29 and 2.31 times larger than the average values measured 

for the background for the beta, gamma and beta-gamma counting systems, 

respectively. However, for the later and longer counts, i.e., 800 minute 

count (Au coin waiting time was T = 7 hrs. + 16 min.) the counting rates 

w 

observed were approximately equal to or less than background and had values 

that were within 3a of the background counting rate, indicative of little 

or no sample contributed radioactivity. The same characteristics is seen 

to hold for the gold foils which had a waiting time T^ = 44 hrs. + 19 minutes 

This suggests that any activity that may hve been present in the gold coin 

or foils, would have had to be quite short lived with respect to the 2.7 

198 

day lifetime of neutron activated Au. 



TABLE I 


SUMMARY OF In FOIL COUNTING RESULTS 


COUNTS 
In it 001 

In it 002 
In it 003 


N B \ N 8y 
O-Om) >gG >BG >BG 


(10m) 'v-BG ^BG ^BG 


(10m) >BG >BG >BG 


SPECTRAL FEATURES 

335, 406, 515, 662 keV (Cs-137) 

340, 505, 662 keV (Cs) 

340, 483, 662 keV 


Combined (40 min.) : 

~BG >BG >BG 


335(1), 414, 662 keV 


Background (40 min.): 


335 (s), 662 keV 
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The ungated gamma spectrum obtained with a 40 minute count of the 

gold coin at 16:45, 12 April 1981, ( > 2% hrs. wait), exhibited no clearly 

distinguishable spectral features not seen in the other ungated gold 

spectra for the longer counting periods and at later times, other than 

137 

the peak at channel 38 (or 340 keV) . The Cs gamma peak and spectral 

features are again sharply defined. 

The spectra obtained during the long term counting measurements for 

137 

the gold coin and foils showed little other than the Cs peak for the 

ungated gamma ray spectra, however the smoothed gated spectra exhibited 

other features, these being the channel 38 peak 340 keV) , as seen earlier, 

137 

and the channels 55-56 peak (y 511 keV) , also seen earlier. No Cs peak 
contribution was observed as can be seen in Figure 9. The same was true 
for the spectra obtained with the multiple gold foil sample, as was shown 
in Figures 10 and 11. Table II summarizes the gold-coin activation measurements. 

198 

The failure to observe the 411 keV Au gamma ray peak in contrast to 
the apparent radioactivation of the Indium foils, resulted in the decision 
to choose indium foils as the detector system for the second series of 
measurements . 

The beta-discriminator was lowered to better accomodate the lower energy 
of the indium beta rays. This led to an approximately doubled beta counting 
rate sensitivity. 

For counting system calibration and set-up, a set of indium foils were 

252 2 2 

activated with both a 5 Ci PuBe source and a small Cf source = 10 n/cm 

seconds) and the samples counted. The system showed excellent sensitivity 

137 

for the indium activity, however, the Cs gamma ray background persisted as 

can be seen in Figures 12 and 13. A careful wash-down with a detergent 

137 

successfully removed most of the Cs activity, especially its effect on 
the gamma ray spectra. The results are shown in Figure 14. 





COUNT RATE PER MINUTE 





FIGURE 11 


MULTI-GOLD FOIL SAMPLE SPECTRUM (GATED) 

(1046, 14 April 1981; to 17 April 1981) 78.9 hrs. 


412 KeV 


I 

I 

i 

) 

i 

i 


511 KeV 


I 

I 

I 


663 KeV 


X - 3-point smoothed data 




19 


TABLE II 


SUMMARY OF GOLD DETECTOR COUNTING RESULTS 


COUNT RATE SPECTRAL FEATURES 




N e 

N y 

N 3Y 


Au 

Coin 40 min. 

2.42BG 

1.29BG 

2.31BG 

40m (ungated) 340 keV, 662 keV 


800 min. 

0.92BG 

1.06BG 

. 78BG 

V36 hrs. (ungated), 662 keV 

Au 

Foils 




340 keV, 511 keV (gated) 


800 min. 

0.86BG 

1.03BG 

• 89BG 

ungated , 662 keV 
gated, 511 keV 
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It should be noted in Figure 14, the gamma spectra of an activated 

In foil and that of background, the peaks observed previously at 662 and 

approximately 340 keV, are absent. This result substantiates that the 

shields were cleaned thoroughly. However, peaks arising from activation 
24 

of Na (15 hrs.) do appear, primarily due to the samples having been 

252 24 

left on the Cf activator set-up for several days, allowing the 4 Na to 

build up further toward its asymptotic limit, while the indium activity 

remained essentially constant. This, it appears likely that the V340 keV 

peak may also be part of the systems response to the background. 

The cleaned and adjusted system was tested, calibrated and then 
transferred to RTS Laboratories for a new series of foil exposures and 
measurements to be undertaken. 

On 13 May, 1981 the discharge was operated for one (1) hour with 

the In foils in place. The In detector foils were then transferred to the 
counter and counted for a 40 minute period. 

The results of the first measurement did not appear to differ signifi- 
cantly from background. Also, during a background count immediately following 
the first measurement, it was noted that the beta and coincidence counters 
would intermittantly count at extremely high rates, resulting in their total 
counts for the period of measurements reaching values of 50 to 100% higher 
than previous values, attempts to identify the sources of the spurious counts 
were unsuccesful at first. The operation of the counting system was 
monitored throughout the day of 14 May and up the point when the high frequency 
for the discharge was turned on. The noise source had apparently disappeared 
and the counting rates due to background was observed to be less than that 
at the regular location of the equipment (the University of Florida). 

In Figure 15 the indium foil spectrum is plotted along with calibration 
22 137 60 

curves and peaks from a Na source (511 keV) , Cs and Co. No obvious 
expression of peaks observed in previous measurements is shown. The information 
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plotted was obtained when the indium foils were used as flux detectors on 
the second round of experiments, 13 May 1981, with the D 2 “discharge (conditions 
were not found to be optimum for the discharge) . 

For comparison, the background spectrum is shown in Figure 16. Many 
of its features appear to have contributed to the other measurements. The 
final figure. Figure 17, shows a peak at approximately 511 keV. This latter 
spectrum was obtained on 15 May 1981, for a set of three In foil detectors, 
which has been exposed in the vicinity of the discharge for a period of 
approximately 1 hour. Even though the total counts, especially the beta 
and coincidence counting rates exhibited a sensitivity to low amplitude, 
high frequency "noise" pulses appearing on the oscilloscope display, the 
gamma ray spectra appear to be relatively free of perturbing aspects, with 
the principal characteristics observed being the peak at 511 keV. 

Conclusions 

The counting rate measurements completed thus far have exhibited certain 
results that have been significantly greater that the average value of the 
observed background counting rates, indicative of radioactivity sources other 
than background. 

Spectral measurements thus far have failed to disclose characteristics 
expected for thermal neutron activated gold, both in terms of peak energies 
and lifetimes. Spectral features characteristic of thermal neutron activated 
indium were observed for some but not all of the indium foil spectra. How- 
ever, peaks have been observed in the spectra at 511 keV, the gamma ray 
energy associated with the positron annihilaiton process. Also, the apparent 
lifetime of the apparent radioactivity appears to be distinctly shorter than 

1QQ 

that due to Au (2.7 days). 

64 

Thermal neutron activated Cu (12.8 hrs.) can yield positrons from 
which annihilation radiation could arise, should copper be a const ituate of 
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the Au coin or foil. However, its lifetime too does not correlate with 
some of the characteristics observed. 

Fast neutron activation through the (n, 2n) reaction can produce 
positron activity in many materials, yet gold’s (n, 2n) reaction products 
are principally electron capture (EC) or beta minus decay species. 

Indium has an (n, 2n) positron emitting product but from its least 
abundant isotope. 

Our efforts this far have posed further questions without providing 
any absolute answers. More carefully controlled measurements are needed 
to assure that the background contributions are better quantitated. A 
series of measurement si, including replicates, need to be performed. 

These should be counting rate, spectral and time dependent. Measurements 
should also be made with the higher resolution germanium detectors. 

The apparent cross sections for the reactions that may have produced 
the observations must be relatively small, causing their measurement to 
be difficult. This is also probably true for the intensity of whatever 
process that may be responsible for the apparent radioactivation. Thus, 
subsequent measurements must be done very carefully in order to optimize the 
exposure times, measurement processes, data storage, retrieval, output in 
graphical and tabular forms. 
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